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ABSTRACT 

The extent to which large scale magnetic fields are susceptible to decay via turbulent 
shredding is important for interpreting the need for in situ large scale dynamos in as- 
trophysics and for observationally inferring field strengths compared to kinetic energy. 
By solving coupled evolution equations for magnetic energy and magnetic helicity in 
a system initiated with isotropic turbulence and an arbitrarily helical large scale field, 
we quantify the decay rate of the latter for a bounded or periodic system. We show 
that energy associated with the non-helical large scale magnetic field decays at least as 
fast as the kinematically estimated turbulent diffusion rate, but the decay rate of the 
large scale helical component depends on whether the ratio of its magnetic energy to 
the turbulent kinetic energy exceeds a critical value given by Mi )C = (fci/fc2) 2 , where 
fci and &2 are the wave numbers of the large and forcing scales. Turbulcntly diffus- 
ing helical fields to small scales while conserving magnetic helicity requires a rapid 
increase in total magnetic energy. As such, only when the helical fields are sub-critical 
can they so diffuse. When super-critical, the large scale helical field decays slowly, at 
a rate determined by microphysical dissipation even in the presence of macroscopic 
turbulence. In effect, the turbulent diffusion of such an initial large scale helical field 
produces small scale helicity whose amplification abates further turbulent diffusion of 
that initial field. Two curious implications are that: (1) Standard arguments support- 
ing the need for in situ large scale dynamos based on the otherwise rapid turbulent 
diffusion of large scale fields require re-thinking since only the non-helical field is so 
diffused in a closed system. Boundary terms could however provide potential pathways 
for rapid change of the large scale helical field. (2) Since Mi jC -C 1 for k\ <C k2, the 
presence of long-lived ordered large scale helical fields as in extragalactic jets does not 
guarantee that the magnetic field dominates the kinetic energy. 

Key words: magnetic fields; galaxies: jets; stars: magnetic field; dynamo; accretion, 
accretion disks; cosmology: miscellaneous 



1 INTRODUCTION 



Many astrophysical sources including galaxies, stars, compact objects, and accretion engines show direct or indirect evidence 
for large scale magnetic fields ( Shukurov|2007 1. The extent to which these large scale fields survive in the presence of in situ 
turbulent diffusion and the conditions that determine their diffusion rates constrains the mechanism of their origin. Do the 
fields result from in situ dynamo generation or could they have been the result of flux freezing from a previous evolutionary 
phase? 

There has been debate over the extent to which 3-D turbulent diffusion of large scale fields is effective and the role 



that the small scale fields play in its potential suppression. The controversy originated in part from 2-D studies ( Cattaneo & 



Vainshtein|1991 1 which seemed to suggest suppression. However, magnetic field lines can interchange in 3-D. This distinction 
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is implicit in the fact that the formalism for computing the isotropic turbulent diffusion coefficient of large scale fields reveals 
a suppression in 2-D that is absent in 3-D ( Pouquet et al.|1976| |Gruzinov fc Diamond|1996 1. The turbulent diffusion of large 
scale magnetic fields has subsequently been studied in terms of an effective turbulent diffusion coefficient for the large scale 
field, scaling this coefficient in terms of some power of the magnetic Reynolds number Rm (Blackman & Brandenburg 2002 
Brandenburg & Subramanian 2005a Brandenburg et al.|2008 1. 



Most work on the diffusion of large scale fields has not distinguished between the diffusion of helical vs. non-helical large 
scale fields. An exception is Yousef et al. (20031, which found that that fully helical fields decay more slowly than large scale 
fields in numerical simulations and discuss this in the context of magnetic helicity evolution. This stimulates a quantitative 
analytic study to understand just how helical the field must be to be resilient to turbulent diffusion. As we will see, the 
suppression of helical field diffusion should be interpreted not as an intrinsic suppression of the turbulent diffusion coefficient 
itself, but as the result of the current helicity correction (Fri sch et al.||1975"| |Pouquet et al.|[f976[ ) to the electromotive force 
which competes with turbulent diffusion. 

One important motivation for this study is the potentially dramatic implications for interpreting the origin of large scale 
magnetic fields in astrophysical rotators (galaxies, disks, stars) as discussed herein. An additional motivation is that jets, 
particular those in active galactic nuclei (AGN) of parsec scale, exhibit Faraday rotation that is consistent with an ordered 

Gabuzda et al.||2008| |201 2\ . This in turn has led some to conclude that the jets 



large scale helical field (Asada et al. 2008 



are necessarily magnetically dominated ( LyutikoT'eT'aL'f 2005 1 . While Poynting flux dominated models of jets ( |Li et al.|200f | 
Lovelace fc Romanova||2003 Lynden-Bell||2006 Contopoulos et al~||2006 1 are plausible, are observed large scale helical fields 



a definitive signature of a magnetically dominated system? 

Although each class of astrophysical source for which the evolution of large scale fields plays a role warrants its own 
focused study, it is fruitful to investigate simplified problems that potentially identify and elucidate basic principles. In this 
spirit, we focus on the specific underlying physics of how long it takes for a large scale magnetic field to diffuse in the presence 
of non-helically forced turbulence when the initial large scale field consists of different helical fractions. We study cases for 
which the initial field strength does not exceed the kinetic energy as the interiors of astrophysical rotators are typically not 
magnetically dominated. 

In section 2 we derive the basic equations to be solved, drawing from previous work on 21st century dynamo theory and 
simplifying the set of equations appropriate for the present problem. In section 3 we solve these equations. We discuss the 
astrophysical implications in section 4 and conclude in section 5. 



2 MEAN FIELD DECAY FROM TURBULENT FORCING IN A CLOSED OR PERIODIC SYSTEM 
2.1 Derivation of basic equations 

Here we derive a system of three differential equations needed to study the decay of large scale magnetic fields of arbitrary 
helical fraction in a closed or periodic box. These are the equations for the time evolution of (i) large scale magnetic helicity, 
(ii) small scale magnetic helicity, and (iii) large scale magnetic energy. From these, we will also construct an equation for the 
evolution of the non-helical large scale field. 



To derive the large and small scale magnetic helicity evolution equations we follow standard approaches ( Brandenburg & 
Subramanian| 2005a Blackm an|2007 1 and start with the electric field 



E = -V$--9 t A, (1) 

c 

where $ and A are the scalar and vector potentials. Taking the average (spatial, temporal, or ensemble), and denoting 
averaged values by the overbar, we have 

E = -V$-^dtA (2) 
Subtracting |2| from |l]) gives the equation for the fluctuating electric field 

e = -V<f> - -<9 t a, (3) 



where <j> and a are the fluctuating scalar and vector potentials. Using B ■ dt A = 9t(A ■ B) + cE • B — cV ■ (A x E), where the 
latter two terms result from Maxwell's equation <9 f B = — cVxE, and the identity A-VxE = E- B — V - (A x E), we take 
the dot product of |2| with B to obtain the evolution of the magnetic helicity density associated with the mean fields 

dt(A ■ B) = — 2cE • B — V • (c<5 B + cEx A). (4) 

Similarly, by dotting Q with b the evolution of the mean helicity density associated with fluctuating fields is 

dt& ■ b = -2ce ■ b - V • (c<j>b + ce^Ta). (5) 
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To eliminate the electric fields from Q and (JsJ we use Ohm's law with a resistive term to obtain 
E = -V x B/c + r/J, 

where J = ^VxB is the current density and rj is the resistivity. Taking the average gives 
E = -e/c - V x B/c + rjj, 



(6) 
(7) 



where £ = v x b is the turbulent electromotive force. Subtracting |7| from (|6| gives 

e = (£ - v x b - v x B - V x b)/c + rjj. (8) 

Plugging |7| into Q and (|8| into |5| and globally averaging (indicated by brackets) to ignore divergence terms gives, for 
the small and large scale contributions respectively 
1 



-&(a- b) = -{£ ■ B) - u M {h ■ Vxb), 
where vm — (r\c 2 /47r). and 



(£ ■ B) - v M {B ■ VxB). 



(9) 



(10) 



To obtain an expression for £ , we use the 'tau' or 'minimal tau' closure approach for incompressible MHD (Kleeorin et al. 



1990} [Blackman & Field 20021; |Brandenburg fc Subra manian 2005a). This means replacing triple correlations by a damping 
term on the grounds that the EMF £ should decay in the absence of B. This gives 



d t £ = (d t v x b) + (v x dth) = "B - ^VxB - £/f, 

T T 



(11) 



where f is a damping time and 
b - Vxb) 



T 

a = — 
3 



4irp 



v ■ V x v) and f3 



— {V >. 

3 X 1 



The time evolution of £ can be retained as a separate equation to couple into the theory and solve, but simulations 
of magnetic field evolution in forced isotropic helical turbulence reveal that a good match to the large scale magnetic field 
evolution in simulations can be achieved even when the left side of (111 is ignored and f = v ^ k , the eddy turnover time 
associated with the forcing scale (Field & Blackmail 2002; Brandenburg & Subramanian 2005b I. We adopt that approximation 
here. Rearranging (111 then gives 



£ = aB- /3VxB, 
Eqs. Q and Q then become 

1. 
2' 
and 

idt(A-B) = a>(B~ 2 ) — /3(B ■ VxB) — vm(B ■ VxB). 



-dt{a-h) = -a(B ) +/3(B • VxB) - u M (b- Vxb) 



(12) 



(13) 



(14) 



Note that the energy associated with the small scale magnetic field does not enter £ above. Therefore it does not couple 
into equations (JoJ) and ( 10 1. It appears only as a higher order hyperdiffusion correction dSubrama nian 2003 1 which we neglect 

h 2 fc 2 

because its ratio to the p term in the EMF is A ° pv -z << 1 . However, upon plugging (12 1 into those equations, the energy 

— 2 2 

associated with the large scale field B does enter. Therefore we need a separate equation for the energy associated with the 
energy of the mean field. To obtain this equation we dot 9 t B = — cVxE with B and ignore the flux terms to obtain 



\d t (B 2 ) 



-c(B ■ VxE) = -c(E • VxB) = {£ ■ VxB) - i/ M ((VxB) 2 ) = a(B ■ VxB) - /3((VxB) 2 ) - i/ M ((VxB) 2 ), (15) 



where the latter two similarities follow from using |7| and ( 12 l and V = 0. 

Eqs. (13 1, (141, and (151 form a set that can be solved in a two scale model as long as {v 2 } is supplied by steady forcing 
such that dt{v 2 } ~ 0, and (v ■ Vxv) remains small. The implications and justification of this latter assumption for present 
purposes will be discussed in more detail after the results of solving the above equations are presented. 



2.2 Two Scale Model and Dimensionless Equations 



To extract the essential implications of the coupled Eqs. (13 1, (14 1, and (15 1 for a closed or periodic system, we adopt 
a standard two-scale model (Blackman & Field 2002 Brandenburg & Subramanian 2005a I and indicate large scale mean 
quantities with subscript "1" and fluctuating quantities with subscript "2." We write the wave number ki > to be that 
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associated with the variation of large scale quantities and that the wave number k 2 » ki to be that associated with the 
variation of small scale quantities. We also assume that k 2 = kf where kf is the forcing wave number at which v 2 — (v 2 ) is 
maintained to be a constant. We also assume that the turbulence is non-helical (i.e. initially driven with ((v ■ Vxv}) = and 
subsequently |(v ■ Vxv}| << |(b • V x b) | / (47rp) , an assumption to be discussed further in section 3.5. Then Eqs (131, (141, 
and ( 15 I then become 

d t Hi -- 



k\H 2 Bl 



9tH a 



and 



dtBl = 



k'lH 2 



4-irp 

4lTp 



2r 2 , 2 rr 
. 2f 2,2rr 



2v M k x Hi, 



2f M k 2 H 2 , 



k%H 2 k\Hi 

4:TYp 



where Hi = (A ■ B) and H 2 = (a • b), B\ 



<B 2 ). 



We non-dimensionlise these equations by scaling lengths in units of k 2 1 , and time in units of r = (k 2 v 2 ) 
the latter equality follows since kf — k 2 in our two-scale approach. We define 



(16) 
(17) 

(18) 
where 



hi 



k 2 Hi 
Anpv 2 ' 



k 2 H 2 
Airpv 2 ' 



v 2 
vuk 2 



and Mi 



Aixpv\ ' 



Eqs. (161, (17 1, and (181 can then be respectively written as 



d T hi 



\h 2 M x 



hi - 



hi, 



(19) 



d T h 2 = 
and 

d T Mi = 



-2 2 
^h 2 Mi + - 



-hih 2 



hi 



-Mi 



2 
RYi 



2 
Rm 



Mi 



(20) 



(21) 



Note that for a fully helical large scale field fei|{A-B)| = (B ). If we choose a positive large scale helicity hi > 0, we can 
drop the absolute value and divide the large scale magnetic energy into a helical fraction fi = (fci(A • B))/((B )) and a non 
helical fraction (1 — /i). Multiplying (19l by ki/k 2 and subtracting it from (21 1, the evolution equation for the non-helical 

Mi — kihi/k 2 becomes 

2 



magnetic energy of the large scale Mi >n h 
2 



d T Mi, nh = -»Mi ltlh 



k 2 k\ 



2 
Rm 



Mi >nh 



(22) 



which has all decay terms and implies a decay rate even faster than that given by the turbulent diffusivity alone when 
hi,h 2 > 0. This will be important in understanding the solution plots that follow in the next section. 



3 DISCUSSION OF SOLUTIONS 



Here we discuss the solutions of Eqs. 1 19 1, ( 20 1, and ( 21 \ for several different cases, assuming that the kinetic energy per unit 
mass is kept steady, and driven by non-helical forcing at kf = k 2 . We identify and derive a minimum helical magnetic energy, 
in units of kinetic energy, required for slow decay. 



3.1 Solutions for fixed initial magnetic energy but varying initial magnetic helicity fraction 

Solutions to Eqs. ( |19[ ), ( |20[ ), and ( |21[ ) are shown in Figs, la-d for Rm = 800 and ki = 1 and k 2 — 5. Each curve in each 
panel represents a solution with a different initial helical fraction /i ; o of large scale magnetic energy, with the initial large 
scale magnetic energy set to equipartition with the kinetic energy, i.e. Mi,o = 1- All cases start with h 2 (t — 0) = 0. The 
six curves of progressively increasing dash spacing correspond to fi t o = 0.95,0.7,0.5,0.2,0.04,0.004 respectively, so that the 
solid lined curves correspond to fi.o — 0.95 and the widest spaced dashed curves correspond to /i t o = 0.004. Fig. la shows 
the time evolution of the large scale helical magnetic energy Mi^ = Tp- hi divided by ^r^, where h c = ki/k 2 is the critical 
magnetic helicity derived in Sec. 3.3 below. The large scale non-helical magnetic field energy Mi iTt /, (Fig lb), and the total 
large scale magnetic energy (Fig lc) Mi are normalized to the initial non-helical magnetic energy Mi >n ho = Mi,o — jj^hi and 
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the initial total magnetic energy Mi,o respectively. The evolution of the non- helical magnetic energy follows Eq. (22 1, derived 
by subtracting ki/k 2 times (191 from (211. Fig. Id shows the time evolution of hijhc. 

All curves of Fig lb show that the non-helical field decays rapidly for all values of /i,o. The minimum decay rates occur 
for /12 = 0. When hi > 0, the rate of decay is even faster than turbulent diffusion, as can be seen from equation (221 in which 
the first term on the right side provides enhanced decay for hi > 0. The rapid decay of the non-helical field in all cases also 
implies that during the slow decay regimes of M\ in Fig. lc, the total field is primarily helical. 

The slow decay regimes in Figs la and lc correspond to decay at a microphysical dissipation rate, determined by the 
last term in Eq. ( 19 1 ) . As the plots show, in these regimes hi(t)/hi, c > 1. When instead hi(i)/hi )C < 1, the decay is much 
more rapid and determined by turbulent diffusion — the penultimate term of (19 1. Correspondingly, when the helical field is 
subcritical right from the start, the helical field rapidly decays. 



3.2 Solutions for fixed initial magnetic helicity fraction but varying initial magnetic energy and scale ratio 



For the solutions to Eqs. ( 19 1, ( 20 1, and (21 1 shown in Fig. 2, we used /i,o = 0.999 for all curves and varied the initial magnetic 



energy Mi,o- We used fci = 1 and ki = 5 with Rm — 800. From top to bottom the curves in Fig 2a correspond to 



20,5,1,0.5,0.1,0.01, respectively. For Fig 2b from top to bottom, the curves correspond to 



Mi, 



A/i.o 
(fcl/fc 2 ) 2 



5,2,1,0.5,0.1,0.01 



(fci/fc 2 ) 2 

respectively. The third curve from the top in each panel corresponds to our analytically derived critical value (see next 
subsection) Mi,o = kih c /k2 = (ki/ki) 2 . This curve marks the approximate demarcation line between slow and fast decay 
curves. For curves with Mi^ ^> (fci/fe) 2 decay is slow (resistive), whereas for Mi,o <§C (ki/ki) 2 decay is fast (unfettered 
turbulent diffusion). 

In Fig. 3 we show solutions to Eqs. ( 19 1, ( 20 1, and ( 21 1 for /i,q = 0.999 but for different values of kijki and Rm than the 
values used in Fig 2. For Fig. 3a, ki = 1, k f = 10, R M = 8000, and for Fig. 3b, ki = 1, k f = 20, R M = 8000, and /i, = 0.999. 



In each panel, the curves from top to bottom correspond to initial energies 



Mi,p 

(fcl/fc 2 ) 2 



5,2,1,0.5,0.1,0.01 respectively. The 



third curve from the bottom in each panel again corresponds to our analytically estimated critical curve bounding the slow 
and fast decay regimes. All curves decay more gradually in Fig. 3b than Fig. 3a. (note the difference in scale of the y — axis 
in the two figures) because both the turbulent diffusion and the resistive diffusion terms for the case considered in Fig. 3b, 
are correspondingly reduced by the smaller value of (fci/A/2) 2 . 



3.3 Derivation of the critical helicity h c 



We now derive the critical helicity h c and the associated critical helical magnetic energy by noting that the slow decay regime 



requires the last term on the right of ( 19 \ to be at least comparable to the sum of the first two terms on the right. For 



large Rm, each of those first two terms is separately much larger than the last term. Therefore the first two terms must 
approximately balance. These same terms also appear in the equation for hi, implying that it too evolves slowly (as seen in 
Fig. Id.) in the slow decay regime. Since initially we always consider hi > and hi = 0, emergence to a slow decay regime 
implies a rapid evolutionary phase (with negligible dependence on Rm) where the buildup of /12 leads to an approximate 
balance between these two terms. But if there is not enough hi , then there is not enough supply of magnetic helicity to grow 
the required hi to establish the slow decay regime. 



We can estimate the needed amount of hi that must be grow by balancing first two terms on the right of either ( 19 \ (or 
p5] |) to obtain that 

h 2 Mi ~~ (ki/kifhi. (23) 

Since Mi ~ kihi/k 2 in the slow decay regime, Eq. ( |23[ ) gives the result that hi ~ ki/ki in this regime. The only possible 
source of hi is hi given our initial conditions, therefore the above value of hi gives the minimum required of hi to achieve 
the slow decay decay regime. That is, we must have hi > h c = kijk% for a slow decay regime of the large scale helical field. 
If ever hi « h c , the large scale field will decay rapidly. Identification of this critical helicity h c is the key to explaining the 
curves shown in Figs. 1,2, and 3. 

The critical magnetic energy associated with h c is simply Mi jC = kih c /k2 = (fci/A/2) 2 an< A this can be substantially below 
equipartition when fci/fe <SC 1. The amount of magnetic energy that decays slowly, Mi tS i ow , is then the difference between 
the magnetic energy contained in the helical field Mi, h = /1M1 = kihi/ki and Mi jC . Dividing this difference by the total 
magnetic energy then gives the fraction of energy that will decay slowly, namely 



Mi iS i ow 
Mi 



M hh - M hc 
Mi 



,0 



Max 



kl 1 



.0 



(24) 



Eq. |24j) shows that most of the initial magnetic energy can decay slowly even if the system is not magnetically dominated or 
maximally helical. For example, we used Mi,o = 1 (equipartition between total magnetic and kinetic energy) and ki/kf = 1/5 
for the solution of Fig. 1 so that a fraction — = fi — i of the initial magnetic decays slowly. 
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3.4 Resilience of helical fields to diffusion is not a reduction in the diffusion coefficient 



The resilience of the helical field to turbulent diffusion as shown above, is the result of the current helicity part of the a effect in 
the language of dynamo theory, not an intrinsic change in the diffusion coefficient f3. This is an important distinction because 
the ft term also appears in the non-helical magnetic energy evolution equation where it is not abated by terms involving 
magnetic helicity, but instead can even be enhanced by them (Eq. |22[ ). The point is that the helical and non- helical large 
scale fields obey different equations. Parametrization of /3 in terms of Rm (e.g. |Blackman fc Brandenburg ( 2002 I; Yousef et al. 



(2003 1) can be misleading in this respect, because /3 itself does not change even when the helical field decays at the resistively 
limited rate. 



3.5 Neglect of kinetic helicity evolution for our choice of initial conditions 

In general the a effect is the difference between small scale current and kinetic helicity, but we have not included an equation 
for the time evolution of the kinetic helicity. If the growth rates of these two helicities were equal so as to keep a ~ 0, then 
the large scale helical field would decay as fast as the non-helical field. We now discuss our justification for ignoring kinetic 
helicity evolution for our specific choice of initial conditions. 

|Blackman fc Fi eld (2004J and |Park fc Blackman| ( |2012| ) showed that the small scale kinetic helicity can grow significantly 
when the system is initiated with fully helical small scale fields. But the small scale current helicity can only drive small scale 
kinetic helicity growth and conserve magnetic helicity by inverse cascading, bringing helical magnetic energy up to larger 
scales. In our present case, any initial magnetic helicity is solely on the large scale to begin with and we now argue that the 
kinetic helicity associated with the small scale is not expected to grow significantly. 

The kinetic helicity growth in the two-scale approximation for a closed system is given by ( |Blackman fc Field|2004| ) 

\dtHX = (~ j ■<^xb)(fc 2 -fci)-^ 2 2 <v-^} + — <w-(jxb)> + ~ <vVx(jxb)>, (25) 
1 \47rp J cp cp 

where v is the viscosity and we assume any non-helical forcing function in the velocity equation does not explicitly contribute. 
We expect that the first term on the right can only grow kinetic helicity if the contribution from either v or j to the correlation 
comes from their helical part. So we can look at this term in two ways, either focusing on the velocity contribution or the 
magnetic field contribution. We consider the latter approach first. If we ignore total divergence (surface) terms and spatial 
derivatives of the cross helicity (v ■ b) (given that the latter evolves only via decay 9t(v • b) = — k\(y + ^m)(v ■ b)), then 



{ijj x b) 9 = (vsdqbs) = (v x }) q and the first term on the right of (25l can be written 2B ■ (v x j)(fc 2 — fci). For a maximally 



helical small scale field j x b = (though Jxb^O^jxB) and so the last two terms of ( 25 1 vanish. If we consider the case 



j • b > 0, then 2B ■ (v x h)k 2 (k 2 - fci) = 2£ ■ Bk 2 (k 2 - fci) and HX could grow positive as quickly as the current helicity k 2 H 2 
since the latter grows at a rate determined by multiplying the first term on the right of Eqn. (Ji)]) by k 2 . 

But since we expect only the helical fraction of the small scale magnetic field to grow kinetic helicity, for non-maximally 

(( ' b\ 2 \ 1/2 
j2 b 2 J < 1 is the helical magnetic 

field fraction at the small scale. The solutions shown in Fig. 1 indicate that k 2 H 2 < kiv 2 /k 2 even when no drain into 
HX is considered. Since we would expect the non-helical turbulent forcing to produce (b 2 ) ~ (« 2 ), we would then have 
/2 = \k 2 H 2 \/ (b 2 ) < ki/k 2 for all /i < 1, and the growth term of HX on the right of (25 I would be less than k\/k 2 times 



that on the right of Q. In addition, when j x b / 0, the triple correlation terms on the right of (251 survive. Since these 
depend only on the non-helical field, we expect them to be decay terms. In addition, an aspect of the kinetic helicity evolution 
that is not well captured in a two-scale theory is that for typical inertial range spectral power laws, the microphysical viscous 
diffusion of kinetic helicity diverges with increasing Reynolds number, unlike that of magnetic helicity (e.g. |Brandenburg fc| 



Subramanian ( 2005a l). This exacerbates the relative importance of microphysical diffusion in the kinetic helicity equation 



compared to that in the magnetic helicity equation. 



Now consider the second possibility that the first term on the right of ( 25 1 survives only when there is a helical velocity 



field with arbitrary magnetic field. The magnitude of the right of (25 1 can then be written 2f v £ ■ Jik 2 (k 2 — fci), roughly a 
factor of the fractional kinetic helicity /„ = j^jj^y slower than time evolution of the current helicity, k 2 H 2 . (The latter again 
evolving at a rate determined by multiplying the first term on the right of Eqn. ([9]) by k 2 .) Since the kinetic energy is forced 
non- helically, f v << I initially and it is likely that £ ■ B would already saturate before significant HX could grow. 

For the above reasons, we therefore expect that HX would not grow significantly to affect our solutions to the specific 
initial value problem presented in the previous subsection. 



4 ASTROPHYSICAL IMPLICATIONS 

Our results show that, when subjected to steady non-helical turbulent forcing of kinetic energy at wavenumber k 2 , the non- 
helical large scale field decays at least as fast as the unfettered turbulent diffusion rate, but the large scale helical field rapidly 
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decays only when the ratio of its energy to the turbulent kinetic energy drops below the critical value Mi, c = (fci/A^) 2 . Above 
this value, the helical field decays on a microphysical resistive time scale. With the caveat that we have not included boundary 
terms or buoyancy, this has several provocative implications. 



4.1 Presence of observed large scale helical fields does not guarantee a magnetically dominated plasma 



Observations of large scale helical fields, such as those detected by Faraday rotation in extragalactic jets ( Asada et aL 2008 



Gabuzda et al. 20121, or inferred in gamma-ray bursts (Yonetuko 20121, are sometimes interpreted to imply that the field 



is force-free and therefore dominates the kinetic energy of the system. Our calculations show that this is not necessarily the 
case: The fact that M\. c <C 1 for k\ <C &2, implies that even significantly sub-equipartition helical fields decay on resistive 
time scales which are typically much longer than dynamical jet time scales. If a jet contained isotropic or quasi-isotropic 



MHD turbulence, perhaps supplied via an instability at the radial interface between jet and ambient medium (Rossi et al 



2008), then the large scale helical field could survive intact and the system would not necessarily be magnetically dominated. 
Although the helical large scale field could appear force-free in the sense that J x B ~ 0, this does not mean that j x B or 



J x b vanish. The non- vanishing of the latter are essential in the derivation of ( 16 1 and 1 17 1, and particularly the appearance 
of H2 (the driving due to the current helicity) on the right sides. 

We did not include any anisotropic velocity such as shear in our calculations. Nevertheless, our results still demonstrate 
that the basic point that mere observation of a helical large scale field does not prove magnetic energy dominance. 



4.2 Rethinking when in situ dynamos are needed to produce large scale fields 



Taken at face value, the survival of helical fields to turbulent diffusion may reduce the essentiality of in situ dynamos in 
systems if boundary terms are unimportant. 

Consider the case of galaxies: A long standing criticism of relying on primordial or protogalactic fields as the primary 
source of galactic fields has been that the mean field would otherwise rapidly diffuse in the galaxy via supernovae induced 



turbulence if this turbulence were unable to also facilitate competitive exponential growth from a large scale dynamo ( Shukurov 



2007). Our calculations provide rejuvenated credence to pre-galactic mechanisms of large scale field production ( |Kulsrud fc 
Zwcibcl 2 008| [Subramanian||2010| |Widrow et al.||2012 I , and specifically those that produce sufficiently strong helical fields 



(Field & Carroll|2000||Copi et al.|2008||Di'az-Gil et al 



2008 Semikoz et al.|2012 1 because only such helical fields would avoid 



diffusing over a galactic lifetime in the absence of boundary terms. 

Most of the energy in large scale galactic magnetic fields resides in non-helical toroidal fields amplified from poloidal fields 
by differential rotation. As long as the turbulent decay time for the non-helical field exceeds the linear shear time, then we 
can expect a predominance of non-helical field in a steady state, even without an in situ dynamo to regenerate the poloidal 
fields. This is because the helical field provides a minimum value below which the toroidal field cannot drop. The toroidal 
field enhancement over the poloidal field would be that which can be linearly amplified in a non-helical field diffusion time. 

One distinguishing signature of a primordial helical field would be that its magnetic helicity would be of one sign and 
would thus not show a reversal across the mid-plane of a rotator. In contrast, an in situ large scale dynamo would be expected 
to produce field whose helicity changes sign across the mid-plane because of the reflection asymmetry of transport coefficients 
that drive the field growth. This leads to predictions for relative signs of large and small scale helicity and their respective sign 
reversals for the sun (Blackman & Brandenburg 20031 which seem to be observed ( |Brandenburg et al.|2011 |, thus providing 
evidence for in situ large scale dynamo action in the sun. We lack such measurements of galactic fields, but the absence of a 
large scale magnetic helicity reversal across the mid-plane would be evidence for primordial galactic fields. 

As mentioned, an important caveat is that our calculations do not include buoyancy or other boundary loss terms that 
could extract large scale helicity at a rate that may still need to be re-supplied from within the rotator. If such terms are 
important, then both helical and non-helical large scale fields would deplete and an in situ dynamo would be needed for 
replenishment. But this shifts the focus from turbulent diffusion in conventional wisdom to that of boundary loss terms in 
assessing the necessity of in situ dynamos. In fact, it may be the boundary loss terms that also facilitate such dynamos in 



the first place by ejecting small scale helicity that would otherwise quench the large scale dynamo (Blackman & Field 2000 
|Vishniac fc Cho|2001| [Shukurov et al.|2006[ |Sur et al.|2007| ). 

Similar considerations regarding the survival of helical fields would apply for the large scale fields of stars and accretion 
disks. A long standing debate over whether the large scale fields that power jets must be produced in situ or survive advection 



in turbulent disks to grow by flux freezing has persisted (Lubow et al. 1994 Lovelace et al. 20091. Our results here would 



suggest that any helical part of these large scale fields would be more resilient to diffusion and more easily advected. Again 
boundary loss terms could change the story in that such terms would both justify the need for a dynamo when it comes 
to ejection of large scale helical fields, while also potentially being essential to its operation by the ejection of small scale 
magnetic helicity. 

Our results would also suggest that, in the absence of boundary terms, fast cycle periods in stars and disks could only 
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depend on a rapid diffusion of non-helical fields, since the helical fields diffuse too slowly. It shou ld be noted that even in a 
closed box, when shear is included, cycle periods can arise with (B } and (A-B) remaining constant (jBlackman fc Brandenburg 



2002 1 . The extent to which such closed volume cycle periods can be fast or remain resistively slow needs more study. 



5 CONCLUSION 



Using an incompressible two-scale theory that couples magnetic helicity evolution and large scale magnetic energy evolution, 
we have quantified the relative decay rates of helical and non-helical large scale magnetic fields subject to isotropic non-helical 
turbulent forcing in a closed (or periodic) system. We identified a critical ratio of the large scale helical magnetic energy to 
the turbulent kinetic energy given by Mi jC = (fci/A/2) 2 above which helical magnetic energy is immune to turbulent diffusion 
and decays at the microphysical resistive rate. In contrast, we find that the non-helical field always decays at least as fast 
as the turbulent diffusion rate and approaches a factor of two faster when the small scale helical field reaches its maximum. 
The calculations herein corroborate the slow decay of helical fields seen in ( Yousef et al.|200"3 1 and provide a more complete 
theoretical explanation by quantifying how much helical field can survive turbulent diffusion and how much decays. 

The physical interpretation of our result emerges from basic principles of magnetic helicity dynamics. For a fixed amount 
of total magnetic helicity, magnetic energy is minimized when magnetic helicity resides on the largest scale. This is the state 
toward which the system would relax in the absence of kinetic forcing (Taylor 1986). If the system is to conserve magnetic 
helicity and the large scale helical field is to turbulently diffuse, then the small scale helical field must gain energy beyond 
than that contained in the initial large scale field. The source of the extra energy is the non-helical kinetic forcing. However 
any growth of small scale magnetic helicity supplies current helicity to the dynamo a coefficient, which in turn regrows large 
scale helical field — the inverse cascade. A quasi-steady state results when there is enough magnetic helicity in the system such 
that the inverse cascade from the small scale sends magnetic helicity back to large scales at a rate competitive with of the 
turbulent diffusion. This state requires an energy in the small scale helical field equal to just fci/fe of the kinetic energy. 

Correspondingly, the initial large scale helical field energy required to supply the necessary magnetic helicity is the further 
reduced fraction (fci/fe) 2 of the kinetic energy. 

All of this leads to a rethinking of the conditions for when in situ dynamos are required in astrophysical objects: If helical 
fields survive turbulent diffusion, then they would only decay on microphysical resistive time scales, obviating the need for 



a source of in situ amplification. This may reinvigorate the potential relevance of primordial helical fields for galaxies ( Field 
fc Carroll|[2000l |Copi et al.||2008| |Dfaz-Gil et al.||2008| |Semikoz et aT]|2012| |Widrow et"aI1|2012| and advected helical fields 



for large scale jets in accretion disks. Lovelace et al. (20091. The fact that large scale helical fields can survive even when 



sub-equipartition with the turbulent kinetic energy also highlights that the mere appearance of a force- free-looking large scale 
field does not prove that the magnetic energy dominates the kinetic energy. 

Our present calculations ignored kinetic helicity evolution, which we argued to be a small contributor to the class of initial 
value problems studied. We also ignored boundary terms and buoyancy, two key ingredients of real systems. Generalizations 
that include these ingredients are of great interest for future work. If, in the context of astrophysical disks and stars, buoyancy 
removes large scale fields without discretion as to their helicity, then in situ dynamos would be needed to sustain large scale 
fields. Work in 21st century large scale dynamo theory has also been evolving toward the perspective that helicity fluxes, if not 
boundary terms, may actually be essential for the operation of large scale dynamos by removing the small scale helicity that 
clogs the evolution of large scale helicity. Thus boundary loss terms for the large scale field may require a dynamo whereas 
the dynamo itself might require the boundary loss terms for the small scale field. 

In this context, note that our initial condition invoked a net helical large scale field without a compensating helical field 
of opposite sign. This circumstance itself could arise from an MHD dynamo only if the compensating magnetic helicity of 
opposite sign were dissipated or lost via a boundary flux. Had we started with equal and opposite helical fields on small and 
large scales subjected to our same non-helical turbulent forcing, the small scale magnetic helicity would inverse cascade and 
annihilate the large scale magnetic helicity rapidly and the total large scale field would decay rapidly. The slow diffusion of 
the large scale helical field that we have studied depends on having a finite net magnetic helicity as a starting point. 
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the six curves of successively increased dash spacing correspond to f\ o = 0.95, 0.7, 5, 0.2, 0.04, 0.004 respectively, (a) Dimensionless large 
scale magnetic helicity, where h c = fci/fc2 (b) Non-helical magnetic energy in units of the initial non- helical magnetic energy; (c) Total 
magnetic energy in units of the initial magnetic energy, (d) Dimensionless small scale magnetic helicity h,2 also normalized to h c . 
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Figure 2. Solutions to Eqs. \Vi\ , l |20| , and for total large scale magnetic energy using a fixed f\fi = 0.999, fci = 1 and = 5 with 
Rm — 800 for various initial magnetic energies Mifi. The third curve from the top in each panel corresponds to our analytically derived 



critical value M 
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Figure 3. Same as Fig 2b. but for &2 = 10 (left panel), &2 = 20 (right panel), and with Rm = 8000 for both panels. 
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